Many-body wavefunctions for normal liquid ^He 
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We present new trial wave-functions which include 3-body correlations into the backflow coordi- 
nates and a 4-body symmetric potential. We show that our wavefunctions lower the energy enough 
to stabilize the ground state energies of normal liquid ^He in the unpolarized state at all pressures 
in agreement with experiment; however, quantitative discrepancies remain. Further, we include 
strong spin coupling into the Fermi liquid by adapting pairing wave functions. We demonstrate a 
new, numerically stable method to evaluate pairing functions which is also useful for Path Integrals 
calculations at low, but non-zero temperatures. 



For the understanding and development of quantum 
many-body systems, liquid helium is a key test case due 
to its experimentally well-characterized low temperature 
phases and the simplicity of the helium interaction. The 
calculation of its thermodynamic and ground state prop- 
erties is one of the most important benchmarks for micro- 
scopic theory and for the quantum Monte Carlo(QMC) 
simulation mcthodsi^Sii^. While QMC can calculate ex- 
act properties for bosonic systems, and has successfully 
been applied to calculate properties for liquid and solid 
''He^, fermion systems are still challenging, since anti- 
symmetry leads to the so-called fermion sign problem. 
In general, no exact method is known which provides a 
solution for an extended fermion quantum system in two 
or three spatial dimensions. 

To overcome the fermion-sign problem, most fermion 
QMC calculations rely on the fixed-node (FN) approx- 
imation where the nodes of a trial wavefunction, ij^T, 
are imposed as a boundary condition on the many-body 
Schrodinger equation with the projector Diffusion Monte 
Carlo (DMC) methodi. Since the nodal surfaces of the 
exact ground state wavefunction are, in general unknown, 
the energies of FN-DMC calculations are higher than the 
exact answer by a small but unknown amount. Progress 
in fermionic Monte Carlo calculations is, therefore, of- 
ten connected to progress in constructing new types of 
trial wavefunctions. Actually, for homogeneous quantum 
systems, only a few types of different trial wavefunctions 
have been successfully applied within QMC, in particular 
for liquid 3HoVd°i". 

Many calculations use the simplest possible nodal 
structure based on a Slater- Jastrow wavefunction 
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where are the positions of the N fermions {i = 
1, . . . N). Antisymmetry is assured by a Slater determi- 
nant, detfei 4>k{^i) of single particle orbitals 4>k{'^), where 
k labels one of the N single particle orbitals and the Jas- 
trow potential Uj[R\ = u{rij), rij = \ri — rj\, takes 
into account correlations {R indicates its dependance on 
all the particle coordinates). Since it is symmetric with 
respect to particle exchange and real, it does not modify 
the nodes of the many-body wavefunction and, there- 
fore, it does not influence the ground state energy within 



a FN-DMC calculation. 

Homogenous quantum systems play an important role 
in understanding many-body correlations, since for inho- 
mogeneous systems, the Hartree-energy typically domi- 
nates the correlation energy. Translational invariance re- 
stricts the single particle orbitals to be plane waves with 
energies less than the Fermi level. This explains why the 
simple Fermi liquid Slater determinant gives accurate re- 
sults, and, at the same time, it explains the difficulty 
in going beyond Slater-Jastrow without violating homo- 
geneity. 

Backflow^- is one possibility to extend this type of 
wavefunction by including many-body correlations in the 
node g-^^d'^d^ using a Slater determinant detki 4>k{tli^) 
where the bare coordinates used in the Slater deter- 
minant of the Slater-Jastrow wavefunctions are replaced 
by dressed "quasi-particle" coordinates 



qf ^ = r. 
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where Yfe[i?] — Vfc vi'^'ij)- additional symmetric 
3-body potential, J/sb, is also added to the Jastrow part 
of the trial function 
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with Wfe[i?] = Vfe X^ij w{rijf^: The resulting backflow- 
3-body (BF-3) wavefunction, 

^BFz = det 0fc(qf ^) exp {-Uj - U^b] , (4) 

ki 

contains the unspecified functions u{r), w(r), and y{r); 
these are determined either by analytical theories or by 
numerically minimizing the expectation value of the en- 
ergy using simple analytical forms with free parame- 
ters. The backflow 3-body wavefunction (BF-3) has pro- 
vided the most accurate results for the electron-gasi2ii2iiS 
and for liquid 3jjq10j11j20j2^ Spin-dependent correlations 
have been considered in RefiSS, but have found to be 
roughly equivalent to backflow. 

However, the calculated liquid '^He energy and mag- 
netic susceptibility is still in disagreement with experi- 
mental estimates^SiSi; the FN-DMC ground state energy 
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is higher by about 260mK. The prediction of the polar- 
ization energy is particularly unsatisfactory: within Vari- 
ational Monte Carlo (VMC), the Slater- Jastrow wave- 
function predicts a polarized fluid at zero pressure^^, in- 
clusion of three-body correlations results in a slight im- 
provement, however, still largely unsatisfactorj»2i,. Back- 
flow mainly affects the unpolarized state and therefore 
improves the estimate of the polarization energy. Only 
within FN-DMC is the unpolarized state stable at zero 
pressure, but the computed susceptibility is still too 
large, and at higher pressure, the unpolarized state be- 
comes unstable.^'''. Because liquid '^He is such an impor- 
tant benchmark system, it is crucial to flnd what corre- 
lations are absent in the backflow-3-body wavefunctions. 

Pairing wavefunctions, where the Slater determinant 
is replaced by the antisymmetrization of a pairing wave- 
function (plji^rj) has been suggested by Bouchaud and 
Lhuillier^^ to overcome this problem. In general, the an- 
tisymmetrization of a pairing functions leads to a pfaf- 
flan, Pfi j(/)(ri, Tj) which reduces to a determinant for a 
M = Q (s-wave) pairing wavefunction. Recently, singlet 
pairing wavefunctions, have been applied to calculations 
on the BEC-BCS crossover of fermionic gases^&SLS^. As 
far as the calculations on '^He are concerned, spin-triplet 
pairing wavefunction with M = 1 (p-wave) are promis- 
ing, but it remains to be shown that they indeed provide 
lower energies than the backflow-3 body wavefunctions. 

It is known experimentally that the ground state of ^He 
is a superfluid with spin triplet pairing. The transition 
occurs around ImK which is an estimate of the energy 
gap, A. We expect that energy differences between the 
true superfluid ground state and the best Fermi liquid 
state are small, of order A ~ ImK and the introduction 
of a superfluid pairing wavefunction alone can not resolve 
the energetic mismatch between theory and experiment 
of order 260mK. Further, one expects that pairing only 
involves states close to the Fermi surface, e.g. states of 
momenta k with kp — 5k < k < kp -\- 5k where k = kp at 
the Fermi surface and 5k ~ 2m,/S./kp defines the typical 
coherence length ^ ~ 5k~^ of the pairing wave function. 
In ^He {kp ~ O.OA^^) one estimates 5k ~ 10"^A~\ and 
the coherence length ^ ^ 100 A exceeds by an order of 
magnitude the size of a typical simulation box, L, as L ^ 
lOA for a system of TV = 66 atoms. Strong correlations in 
helium in changing the bare mass into the effective mass 
hardly modify this conclusion. On the other hand size 
effects (for a non-degenerate ground state) will favor the 
normal state since the energy of a weakly bound state 
is increased when the size of the system is decreased. 
For these reasons, we mainly limit the discussion to the 
normal state. 

In this paper, we propose new types of trial wavefunc- 
tions with explicit many-body character and use them 
to calculate the unpolarized ground state energy of '^He 
at the maximum, intermediate, and the minimum densi- 
ties of liquid ■^He. First, we introduce a three-particle 
backfiow and 4-body correlations (3BF-4), as a natu- 
ral extension of the backflow-3-body (BF-3) Fermi liquid 



P (nm ^) 


wavefunction 






Edmc{K) 


16.35 


BF-3 


-2.201(6) 


23 


-2.417(1) 




3BF-4 


-2.284(3) 


14 


-2.438(1) 




S3BF-4 


-2.294(3) 


14 


-2.432(4) 


19.46 


BF-3 


-1.775(5) 


26 


-2.155(5) 




3BF-4 


-1.905(4) 


20 


-2.174(3) 




S3BF-4 


-1.945(4) 


20 


-2.182(4) 


23.80 


BF-3 


-0.055(1) 


65 


-0.77 (2) 




3BF-4 


-0.272(7) 


53 


-0.834(2) 




S3BF-4 


-0.340(5) 


49 


-0.861(4) 



TABLE I: Variational (£„) and FN-DMC energies (Edmc) in 
K/atom for a system of V = 66 unpolarized '^He atoms with 
periodic boundary conditions at two different densities; is 
the variance per atom; BF-3 are results of the backflow-3body 
wavefunction, 3BF-4 are the results presented here using an 
additional 3-body correlation for the backflow and an addi- 
tional symmetric 4-body term, and S3BF-4 is the pairing form 
described below. 

wavefunction. As in the case of backflow, the Slater- 
determinant is a function of quasi-particle coordinates, 
but a three-body correlation is added to the original 
backflow coordinate 

qr^ = qr + Q.D, (5) 

where Qi[^] is a tensor with generic elements 
Qk[R])ai3 = VfcaVfe/j li^j) + 5ai3^l J^ij I'irij), and 
Dfc[i?] — Vfe^j^ d{rij) is a vector. Furthermore, a sym- 
metric 4-body potential is added to the Jastrow part of 
the trial function: 

C/4B = ^X,PiXi (6) 

i 

where Pi[R] is a tensor and a vector, analogous to 
Q and D. The new trial wavefunction including 3-body 
backflow and a 4-body potential (3BF-4) is: 

V'3BF4 = det M<ir'') exp[-C/j - U^b - U4b]- (7) 

ki 

As in previous studies, we use gaussian functions as a 
basis for the 1-D functions arising in the wavefunction 
(such as q,q' and d), and minimize a combination of en- 
ergy and variance of the trial wavefunction within VMC. 
Results for a system of iV = 66 unpolarized '^He atoms 
interacting with the recent Aziz-SAPT2 potentialS2i2£ 
with periodic boundary conditions are shown in table 
Si. To reduce flnite size effects we have also performed 
twist averaging'^^ and we have corrected size-effects in 
the potential energy. Results for different polarizations 
are given in table II. For comparision we have added the 
twist-averaged results of the backflow 3-body wavefunc- 
tions of Ref.^^ in table II and the corresponding results 
for periodic boundary conditions in table I. 

The three densities 16.35, 19.46, and 23.80 atoms 
nm~^ correspond to pressures equal to 0, 8.4, and 35.3 
bar. The magnetic susceptibility, x/C*? has been esti- 
mated in table III by assuming E = Eq + C^/(2x/C) -I- 



p 


c 


E^{K) Edmc{K) 
BF-3 BF-3 


E^{K) Edmc(K) 
3BF-4 3BF-4 


16.35 




0.242 
0.485 


-2.1633(9) -2.3586(8) 
-2.1614(9) -2.3548(10) 
-2.1753(9) -2.3433(8) 


-2.241(1) -2.3802(4) 
-2.242(1) -2.3743(1) 
-2.248(1) -2.3633(6) 


19.46 




0.242 
0.485 


-1.7469(12) -2.0685(13) 
-1.7499(12) -2.0623(13) 
-1.7984(12) -2.0807(11) 


-1.863(1) -2.0972(9) 
-1.873(1) -2.0948(6) 
-1.905(1) -2.1046(8) 


23.80 




0.242 


-0.0184(20) -0.6612(20) 
-0.0229(21) -0.6619(21) 


-0.203(2) -0.7094(9) 
-0.222(2) -0.7093(11) 



TABLE II: VMC and DMC energies of the BF-3 wavefunction 
of Ref. ^ and of the 3BF-4 wavefunction using twist averaged 
boundary conditions. Notations and units as in table I; C, is 
the polarization. 
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(nm-^) 


BF-3 


X/C{K-^) 
3BF-4 


pcxp 

(nm-^) 


Xcxp/C 


16.35 


8(2) 


4.5(5) 


16.37 


3.0 


19.46 


5(2) 


8(3) 


19.44 


4.0 


23.80 




|X/C|>14 


23.45 


6.1 



TABLE III: Magnetic susceptibility x/C in where C is 

the molar Curie constant, estimated from the DMC calcula- 
tions using the 3BF-4 wavefunction; experimental values— at 
densities p^xp are given for comparision. 

a4("* where C is the spin polarization. At the equilib- 
rium density, inclusion of 3 and 4-body correlations re- 
duce the magnetic susceptibility in agreement with the 
experimental valueiS. Away from equilibrium density, 04 
becomes more important, and the estimate of x/C' at the 
intermediate density is less reliable. Close to crystalliza- 
tion, the energy variation between the unpolarized and 
the partially polarized state for the 3BF-4 wavefunction 
is of the order of the uncertainty of the calculation. We 
have not performed calculations for the 100% polarized 
phase since backflow is already much less important than 
in the unpolarized phaseSS; similarly 3-body correlations 
mainly affect the unpolarized state24. Therefore, we do 
not expect significant corrections to the energies of the 
complete polarized state of the BF-3 trial function. 

We have tried also different tensorial forms for a 4- 
body trial function, a 5-body trial function, as well as 
4-body and k-dependent corrections to the quasi-particle 
coordinates without significant improvements (< 5mK). 
We cannot exclude the possibility that more effective op- 
timization with more flexible functions would lead to sig- 
nificantly lower energy. This possibility, together with 
the application of the released node technique'^^ is under 
current investigation. 

We have also considered pairing wavefunctions. Sim- 
ple functional forms (Gaussians) for the pairing functions 
(plus the Jastrow potential) have led to energies consid- 
erably higher than those from the BF-3 trial functions. 
Shell effects lead to large size effects and make a direct 
comparison of energies of triplet and singlet pairing wave- 
functions problematic without a careful extrapolation to 
the thermodynamic limit. Since triplet pairing functions 
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did not show significant improvements over the singlet 
pairing, we restrict the pairing functions to spin-singlet 
pairing: the pfaffian reduces to a determinant and the 
antisymmetric part of the wavefunction is: 

det0s(r»T - Tji) (8) 

where r^j {^ii) represents the coordinates of particle i 
with spin-up (down). Using a Gaussian form for (t)s{r), 
we have that the minimum energy orbital has a width as 
large as the simulation cell, L. In this case, the orbitals 
have to be periodized, which is most easily done using a 
Fourier sum. Optimization of the Fourier coefficients led 
to the limiting form of the Fermi liquid: 

k<kF 

We conclude that pairing of bounded singlet /triplet pairs 
does not improve the energy significantly, consistent with 
the small pairing energy in superfiuid '^He.'^^ 

To go beyond the Fermi liquid trial function, we can 
use the singlet pairing wavefunction as a way of intro- 
ducing strong spin-correlations: 

(j)^(r) = J2 e'''"'[i+''('''le"''('^) (10) 

in the determinant of Eq. ijHJi. Here 7y(r) and i'{r) are 
localized functions vanishing at large particle separation. 
They introduce modulations and correlations in the pair- 
ing of opposite-spin atoms without affecting the long- 
range correlations of the Fermi- liquid determinant. They 
do not describe a bound state, so we do not expect them 
to describe a superfluid state, but they might be a precur- 
sor to superfluidity. However, the correlations between 
unlike spins turn out to be much stronger than introduc- 
ing a spin-dependent backffow or Jastrow potential which 
have not been successful in lowering the energy. We also 
include backffow, 3-body backflow, 3-body, and 4-body 
potentials as discussed above: 

^S3BF4 = det0.(q?f ^-q,\^^)e-^-^--^-. (11) 

As shown in table I, the energies are lowered by this wave- 
function. By construction, this wavefunction does not 
affect the energy of the completely polarized state and, 
therefore, it stabilizes the unpolarized phase. We have 
not performed twist-averaging, but it is likely that the 
relative gain of the energy compared to the 3-body back- 
flow wavefunction survives in the thermodynamic limit, 
but a more detailed study of the finite size extrapolation 
of this wavefunction remains to be done. In Ref."^, it was 
estimated that the 3-body potential together with the ef- 
fect of higher order many body terms in the interparticle 
interactions raise the energy by ^ 140mK, and have to be 
added to the computed results when compared to the ex- 
perimental energySS, —2A81K at p = 16.35nm~'^. Even 
if the absolute energies are therefore still in disagreement 
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with experiments by 240mK the agreement with the 
compressibiUty is considerably improved. Further, we ex- 
pect a reduced magnetic susceptibility compared to the 
3BF-4 wavefunctions, closer to the experimental values. 
At high density 10 — 20mK might be still missing to sta- 
bilize completely the unpolarized phase. However, size 
effects are of the order of 30mK, and a clear prediction 
would involve an extensive finite size study. 

We briefly explain our method to evaluate the deter- 
minant of the pairing function in Ea. Hllf) . A straight- 
forward evaluation of the determinant of the matrix 
^ij = i'sij'ii fj) is numerically unstable in the limit where 
the wavefunction approaches a Fermi liquid, Eq. (jSJ, be- 
cause it becomes ill-conditioned; a direct evaluation of 
the determinant is dominated by round-off errors as the 
number of particles, N, increases. Let us consider the ex- 
pansion of the pairing functions, into eigen-modes, ifk{r), 

(t>s{n,rj) ^^dkifl{ri)ifk(rj) (12) 

k 

where dk is the occupation number of eigenstate k. For a 
Fermi liquid, the summation is dominated by N modes, 
k = 1, . . . , A''. We introduce the matrices of the Slater- 
determinants of the Fermi Liquid, Lik — fKri), Rki — 
fk{fi), and a diagonal matrix Dkk' = Sk.k'dk, and write 
Ea. ()12|l as a matrix, 

* = LDR + M = LD^/^ (^1 + D^/'^R (13) 
with My = '^lk>N'^k'Pl{ri)^pk{rj), and M = 



L-i£)-i/2jif£)-i/2^-i^ Now, the determinant of 
Ea. p3(l can be written as: 

det$ = detLdetDdeti?det(l + M) (14) 

and involves only the usual Slater determinants, detL, 
det R and their inverses which are well-behaved. The ma- 
trix M is typically much more rapidly decaying in space 
and a numerical evaluation of the terms on the rhs of 
Eq. (|14() is numerically stable, including all limiting cases. 
This method is not restricted to spin singlet / triplet pair- 
ing, but it is also useful to calculate finite temperature 
density matrices well below the Fermi-temperature where 
the same problem of extremely ill conditioned determi- 
nants occurs. 

In conclusion, we have found 3-body and 4-body cor- 
relations improving the energy and the magnetic suscep- 
tibility of liquid "^He. We expect that the new terms will 
be important to calculate the Fermi liquid parameters 
in strongly correlated liquids such as liquid "^He and the 
electron gas, the polarization transition of the electron 
gas close to Wigner crystallization and related systems. 
The method to calculate pairing function close to the 
Fermi liquid should further allow the study the very di- 
lute limit of a BCS gas. 
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